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Abstract 


Emission lines of [0 III] at 52 pm and 88 pm and of [N III] at 57pm in 
the nucleus of the galaxy M82 have been observed from the Kuiper Airborne 
Observatory with the facility's cooled grating spectrometer. The [N III] 
line has not been previously detected in any extragalactic source. The 
fluxes in the lines indicate -4*1 0 7 M© of ionized gas and a large population 
of massive stars (equivalent to 5*10 5 08.5 stars), sufficient to power the 
infrared luminosity of the nucleus. We use the 52 to 88 pm line intensity 
ratio to find an average electron density of 210 ± 75 in the nucleus; this 
is 10 to 100 times lower than values typically observed in individual 
compact HII regions in our Galaxy. The relative line strengths of the 
[0 III] and [N III] lines imply an N ++ /0 ++ ratio of .45 ± .1, significantly 
lower than is measured by the same method in individual HII regions at 
similar galactocentric distances (S400 pc) in our Galaxy. This lower N* + /0 ++ 
ratio may be due to a lower N/0 ratio, higher stellar temperatures, or both, 
in M82. At spectral resolutions of -90 km/s, all three line profiles are 
similarly asymmetric. They can be well fitted by two Gaussian distributions 
with widths of -150 km/s and central velocities of -110 and -295 km/s, 
bracketing the systemic velocity of the nucleus of -210 km/s. Within 
uncertainties, both the N ++ /0 ++ ratio and the electron density are the same 
for both Gaussian components; this indicates no major large-scale gradient 
in either quantity within the nucleus. 



I. Introduction 


The nucleus of the irregular galaxy M82 is thought to be an active site 
of recent star formation. Models based on previous observations predict a 
greater proportion of early-type stars there than in our own Galaxy (Rieke 
et al. 1980; Kronberg, Biermann, and Schwab 1985). Since the far infrared 
[0 III] and [N III] lines we have observed are emitted primarily in regions 
ionized by early-type stars, and since these lines suffer only minimal 
obscuration, their fluxes are a sensitive probe of the presence of 
early-type stars throughout the nucleus. In addition, the ratios of these 
line intensities allow the electron density and N ++ /0 ++ ratio in the 
interstellar medium of M82 to be determined and compared to those in our 
Galaxy. 

Several previous measurements of M82 show similar, highly asymmetric 
line profiles. The CO 1-K) line measured with a 50" beam (Young and Scoville 
1984), and the Ha (Heckathorn 1972) and [Ne II] 12.8 pm (Beck et al. 1978) 
lines in synthesized 1* beams all have widths of -300 km/s (FWHM), mean 
velocities roughly equal to the systemic velocity of the galaxy (-210 km/s), 
and peak intensities at substantially lower velocities (-110 km/s). 

Previous far IR measurements of the 52 and 88 pm [0 III] lines (Duffy et 
al . 1984; Watson et al. 1984; Lugten et al. 1985), the 158 pm [C II] line 
(Crawford et al. 1985), and the [0 I] 63 pm line (Watson et al. 1984) may be 
consistent with the CO, Ha, and [Ne II] lines but have poorer 
signal-to-noise. 
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In section II, we present new measurements of the [0 III] line? at 
higher spectral resolution and improved signal-to-noise which show similar 
profiles to the CO, Ha, and [Ne II] lines. We also report the first 
measurement in an extragalactic source of the [N III] line, which has a 
similar line profile. From these measurements we derive the electron 
density and N ++ /0 ++ ratio in the nucleus. In section III we discuss these 
results, estimate the number and temperature of the exciting stars in the 
nucleus, and discuss physical interpretations of the asymmetric [0 III] and 
[N III] line shapes. Our conclusions are given in section IV. 

II. Observations and Results 

M82 was observed with the facility cooled grating spectrometer (CGS) on 
the 91.4 cm telescope of NASA's Kuiper Airborne Observatory. The CGS has a 
liquid helium cooled echelle grating 45 cm in length which disperses light 
onto a detector array. For these observations, 6 discrete Ge:Ga 
photoconducting detectors were used to measure intensities at adjacent 
wavelengths. Different wavelengths are observed by changing the angle of 
the echelle. Resolving power X/AX varies from -1700 to -6000 depending on 
echelle order and aperture size. The CGS is described in detail by Erickson 
et al. (1984, 1985). 

The 51.815 pm and 88.356 pm transitions of [0 III] were observed on the 
night of 6 February 1985; the 57.330 pm [N III] line was measured 4 February 
1985. The infrared beam was centered at ct^Q = 9 h 51 m 42.8 s ; - 

69°54'59" using a nearby guide star. This position is nearly half way 



between the two 10 pm continuum peaks of Rieke et al. (1980). Estimated 
pointing error was ± 5 arcsec. The chopping secondary mirror was driven 
perpendicular to the major axis of the galaxy, with an amplitude of 4'. Our 
effective beam diameter of 48" included the roughly 30" by 15" region of 
ionized gas in the nucleus (Beck et al. 1978). 

M82 was observed at spectral resolving powers of 2900 (88 ym) , 2600 
(57 ym), and 3300 (52 ym). At these resolving powers, the widths of the 
Doppler-broadened spectral lines exceed the total bandpass of the 
spectrometer’s 6 detectors. Therefore, 3 or 4 different echelle positions 
were used to provide contiguous wavelength coverage over each line and its 
nearby continuum. 

Raw spectra were corrected for absorption by atmospheric water vapor 
and instrument response by dividing them by calibration spectra of the 
Kleinmann-Low (KL) nebula taken the same night at the same or a nearby 
wavelength. The absolute flux calibration was then obtained by multiplying 
these ratios by the KL flux as measured by Erickson et al. (1981). The 
continuum of KL was observed at only one of the echelle positions at which 
each line in M82 was observed. Calibration spectra at the other (nearby) 
wavelengths were generated by correcting the measured KL spectrum for 
changes in the wavelength-dependence of the atmospheric transmission and the 
instrument response. The water vapor correction was calculated assuming 7 
precipitable microns of water in the line of sight as inferred from fitting 
a narrow unresolved atmospheric water line measured on each flight with the 
CGS. Errors introduced by uncertainties in this correction are less than a 
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few percent. The relative instrument response at different wavelengths was 
measured on the ground immediately after each flight using a blackbody 
calibration source. These corrections are 15? or less at 88 and 57 pm, but 
as large as a factor of 1 .75 at 52 pm due to variations in the transmission 
of the instrument caused by interference effects in the CGS long-pass 
filter. Calculated diffraction losses are less than 10? and somewhat 
uncertain, so we did not correct for them. The final spectra are shown in 
figures 1 - 3. 

Both [0 III] lines and the [N III] line have asymmetric profiles which 
are qualitatively quite similar. To obtain flux estimates, we fitted each 
spectrum with a curve consisting of a constant continuum plus two superposed 
Gaussians. The free parameters in each fitted curve were the continuum flux 
density and the intensity, half-width and peak wavelength of each Gaussian. 
The fits are good for all three lines (x 2 per degree of freedom £1) and are 
also shown in figures 1 - 3- 

In table 1 we list velocities and fluxes for each line as a whole and 
separately for each Gaussian component in each line; we refer to the two 
components as A and B. Radial velocities were calculated from the centroid 
of each Gaussian, using rest wavelengths from Moorwood et al. (1980). 

Quoted velocity uncertainties are largely due to systematic uncertainties 
(estimated at ±25 km/s) in the wavelength calibration, which is based on 
laboratory observations of water lines and in-flight observations of the 
Trapezium. For the 57 pm line, the uncertainty in the rest wavelength is 
also significant (±17 kra/s). These systematic uncertainties do not affect 
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the velocity separation of components A and B, which is '185 km/s in all 
three lines (table 1). Quoted uncertainties in the line component widths 
represent the statistical uncertainty in fitting the curves to the data. 

Flux uncertainties in table 1 are statistical only; systematic 
uncertainties are estimated at an additional ± 20$ of the fluxes. The ratio 
of fluxes in components A and B, which is relatively unaffected by 
systematic uncertainties, also does not differ significantly from line to 
line. Absolute line fluxes have been corrected for foreground extinction in 
M82 by extrapolating the value of t - .30 at 33 pm (Houck et al. 1 98 4 ) using 
a 1/A wavelength dependence. These corrections are 21, 19, and 1 2% at 52, 
57, and 88 pm respectively; the line flux ratios are thus not strongly 
affected by extinction corrections. 

Our observed continuum flux densities are 20 % to 4 0 % greater than the 
broad-band photometric measurements of Telesco and Harper (1980), who used a 
50" beam. Our 52 pm and 88 pm line fluxes differ by 15$ and 35$ 
respectively from our previous measurements made with a 37" beam (Duffy et 
al . 1984). These differences are all within estimated combined 
uncertainties. However, our 88 pm line flux is roughly half the value 
reported by Watson et al. (1984) with a 44" beam. This discrepancy is not 
understood. 

The relative line strengths of the 52 pm and 88 pm lines can be used to 
determine the average electron density in the emitting gas (0 ++ regions). 

We assume the transitions are excited by collisions with electrons and are 
optically thin. Using a 5 - level model 0 ++ atom (Rubin 1984) with atomic 
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constants from Mendoza (1983), we find n . 210 ± 75 in the nucleus. 

Assuming that some of the systematic flux uncertainties (e.g. that due to 
error in the assumed flux of the calibration source) cancel when taking line 
flux ratios, we used the statistical plus half the systematic uncertainty as 
the total uncertainty when calculating the density. Since the spectra can 
no doubt be well fitted by mathematical functions other than the two 
Gaussians we have used, we do not claim that these Gaussians represent 
physical sources. However, because the relative flux in the two components 
is indistinguishable for the two [0 III] lines, our spectra show no evidence 
for a large-scale gradient in electron density within the nucleus. 

At the low density seen in M82, the [0 III] line ratio is more density 
sensitive (Saraph and Seaton 1970) than the [S III] 33 • ^ ym/18.7 vim line 
ratio measured by Houck et al. (1984). They measured an average electron 
density in S ++ regions of 120*^80 cm”®* which is lower than but consistent 
with our value. Their observations do not have sufficient spectral 
resolution to detect asymmetries in the line profiles. 

We calculate N ++ /0 ++ abundance ratios from the relative strengths of 
the [N III] and two [0 III] lines and the atomic models of Rubin (1984). We 
find a ratio of .45 ± .1 for the whole nucleus, again including estimated 
systematic errors as discussed above. Like the electron density, the N ++ /0 ++ 
ratio does not show large variations across the nucleus. 

III. Discussion 
1 . Properties of Exciting Stars 
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Models of M82 based on a wide range of observed properties indicate 
that it has a large population of early-type stars (Rieke et al. 1980; 
Kronberg, Biermann, and Schwab 1985). We use observed infrared line fluxes 
and H II region models to independently estimate the average effective 
temperature, number, and total bolometric luminosity of these stars. 

The H II region models of Rubin (1985) show that the luminosity in the 
[0 III] 52 pm line relative to that in Bra depends strongly on the effective 
temperature of the exciting stars but only weakly on the interstellar 
electron density or the number of ionizing photons. Using the observed Bra 
flux of the nucleus Willner et al. (1977) and our 52 pm line flux, Rubin's 
models imply a stellar temperature of 35500 ± 300 K (spectral type 08.5) for 
uniform density H II regions having 0/H ratios within a factor of 2 of solar 
( 6 . 76 . 10 ~ 4 ). 

Our observed [0 III] line luminosities and those predicted for model H II 
regions having this stellar temperature imply that the [0 III] lines are 
excited by the equivalent of 5 - 1 0 5 08.5 stars in the nucleus of M82. By 
comparison, Mezger (1982) estimates that 0 stars in the entire Milky Way 
produce 3*1 0 s * Lyman continuum photons per second, equivalent to only 2*1 0 s 
08.5 stars (Panagia 1973). 

The early-type stars needed to produce the observed [0 III] lines also 
account for the observed bolometric luminosity ( 4 - 1 0 1 0 L©; Telesco and 
Harper 1980) of the nucleus. At 5.4 -10** L© each, the 5 * 1 0 s 08.5 stars 
required to produce the [0 III] lines have a total bolometric luminosity of 
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2.5-10 10 L©, nearly the observed value. Most of this luminosity appears in 
the far IR because the short wavelength starlight is absorbed and reradiated 
by dust. 

2. Electron Densities 


The average electron density we measure in M82 (210 ± 75 cra~3) i S much 
less than densities measured in individual Galactic H II regions by the same 
method; those values are typically in the range 1,000 to 20,000 cm“3 (Lester 
et al. 1983). However, the higher densities found in our Galaxy may result 
primarily from an observational selection effect. Radio observations 
(Mezger and Smith 1975; Smith, Biermann, and Mezger 1978) indicate that only 
10 - 20$ of the 0 stars in our Galaxy are surrounded by compact (n e >-io* cm~^) 
H II regions; the remaining 80 - 90$ have "extended low density" (ELD) H II 
regions with typical densities of -10 cm~3 (Mezger 1978). The Galactic H II 
regions most often observed, those with high surface brightnesses and large 
densities, are thus atypical; large-beam observations of our Galaxy would 
measure a lower electron density than is seen in individual compact H II 
regions. 

Houck et al. (1984) present a simple model of the interstellar medium 
in M82 which corrects for this selection effect and accounts for the 
observed electron density. They conclude that there is proportionately less 
high density ionized gas in M82. Our observed electron density for M82, 
which agrees with that of Houck et al . , supports this conclusion. This lack 
of high density gas may be a result of the high supernova rate (0.2 per 
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year) in M82 (Kronberg and Sramek 1985); winds from these supernovae may 
blow significant amounts of interstellar gas out of the nucleus of M82 and 
thus reduce the gas density (Chevalier and Clegg 1985; Wyse and Silk 1985). 

3. Mass of Ionized Gas 


Our measurement of the electron density in M82 allows us to improve the 
estimate of Willner et al. (1977) of the mass of ionized gas implied by 
their Bra line flux. To first order, the Bra flux per unit mass of ionized 

gas lj Bra /M HII in H 11 regions is proportional to the interstellar electron 

density and independent of the effective temperature of the exciting 

star(s). For an average electron density of 210 cm~3 > Rubin's H II region 

models and the observed dereddened Bra luminosity (8.5 *10 6 L©; Willner et 

al. ) give an ionized gas mass of 4*10 7 M© in the inner 30" of M82. This is 

naif Winner's value of 7*10 7 Me, which was obtained assuming a volume 

filling factor for ionized gas f = (n ave /n rra3 ) 1 /2 =1. 

e e 

4. N/0 Ratio in M82 


The average N ++ /0 ++ ratio (.45 ± .1) we measure in M82 can be 
compared to those obtained by the same method for individual Galactic H II 
regions by Dinerstein et al. 1984. They find values from 0.2 to 1.2 with a 
tendency to increase towards the Galactic center. Our N ++ /0 ++ ratio 
pertains to the central 400 pc of M82; the two measured H II regions within 
this distance from the center of our Galaxy, Sgr A and G0.5-.1, have higher 
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N ++/0++ rat i oa (0 .7^*2 and 1.1 ± .2, respectively). These Galactic values 
would be even higher (by -50?) if derived using the up to date Rubin model 
atom used here. 

0 

The lower N ++ /0 ++ ratio in M82 may reflect a lower N/O ratio. Low N/0 
in M82 is expected from the models of Rieke et al. (1980) and Kronberg, 
Biermann, and Schwab (1985) which indicate a large population of high-mass 
(M/M© > 20), hot (T > 35000 K) stars and an extreme shortage of stars below 
5-10 M®. Numerous high-mass stars imply a large oxygen abundance, since 
supernovae from these stars are the principal source of oxygen in the 
interstellar medium (Weaver et al. 1978). Few lower-mass stars may 
additionally imply a lower nitrogen abundance, since the main source of this 
element is stars of 5 - 10 M© (Wyse and Silk 1985). 

However, another effect can influence the observed N ++ /0 ++ ratio. The 
models of Rubin (1985) indicate that although 0 and N are both almost 
completely doubly-ionized in H II regions excited by very hot stars, a 
larger fraction of N than of 0 is doubly-ionized by cooler stars (figure 4). 
This means that the N ++ /0 ++ ratio approaches the actual N/0 ratio in H II 
regions excited by hot stars, but is significantly greater in 
lower-excitation regions. The size of this effect is largely independent of 
density and clumping of the ionized gas, and depends weakly on the 
interstellar 0/H ratio (figure 4). Thus the lower N ++ /0 ++ ratio in M82, 
while suggestive of a lower N/0 ratio, may instead be completely due to a 
predominance of stars hotter than those exciting Sgr A or G0.5-.1. 

Because N ++ /0 ++ ratios depend very strongly on stellar temperature, and 


- 12 - 



because these temperatures are not well known in either M82 or the Galactic 
center, our N ++ /0 ++ ratios and those of Dinerstein et al. are only upper 
limits to the corresponding N/0 ratios. These limits are consistent with 
each other and with the N/0 ratio of 0.2 found for the Orion region from 
models based on IR line measurements (Simpson et al. 1986). 

5. Sources of Line Emission in M82 


It is clear from the measured line profiles that neither 0 ++ nor N ++ 
emission is uniformly distributed throughout the nucleus. Two possible 
explanations are (1) the stars which ionize this gas (those having effective 
temperatures above -30000 K) may be nonuniformly distributed, for example in 
a ring in the galactic plane (which we see edge-on) or (2) the gas and dust 
densities are nonuniform and peak away from the center of the nucleus, 
perhaps because of strong stellar winds flowing out from the center. Either 
possibility can account for the two maxima in the 10 pm emission (Rieke et 
al . 1980) by producing two maxima in the temperature or the column density, 
respectively, of the emitting dust. The location of the 10pm peaks and the 
Beck et al. (1978) rotation curve give velocities for these peaks of 130 and 
250 km/s, which would be consistent with those of the Gaussians fitted to 
the [0 III] and [N III] lines (about 110 and 295 km/s) if the coordinate 
scale of the 10 pm map were shifted 1 - 2" relative to that of the rotation 
curve. This agreement suggests that the two components fitted to the 
observed line profiles could correspond to the two 10 pm emission peaks. 
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IV. Conclusions 


1) Fluxes in the [0 III] 52 and 88 pm lines indicate excitation of the 
interstellar medium in the nucleus of M82 by the equivalent of 5*10 s 08.5 
stars. The bolometric luminosity of these stars roughly equals that of the 
nucleus, suggesting that they are its primary energy source. This is 
consistent with the picture of M82 as an active site of recent star 
formation. 

2) The electron density of 210 cm - 3 derived from the 52 to 88 pm line 
ratio is consistent with the value obtained by Houck et al. (1984) from 
lower resolution measurements of the [Sill] 33.4 and 18.7 ym lines. These 
densities are much lower than those measured by the same method in 
individual H II regions in our Galaxy. We support the conclusion of Houck 
et al . that there is proportionately less high-density gas in M82 than in 
our Galaxy. 

3) Our determination of the electron density in M82 allows us to 
recalculate the mass of ionized gas in the nucleus implied by the Bra flux 
of Willner et al. (1977). We find 4«10 7 M© of ionized gas, roughly half 
Winner's estimate. 

4) The N ++ /0 ++ ratio derived for the nuclear gas from our ionized line 
fluxes is less than those derived in the same way for two H II regions at 
comparable galactocentric radii in our Galaxy. This may result from (1) a 
lower N/0 ratio in M82, possibly due to an overabundance of oxygen caused by 
frequent supernovae and/or an underabundance of nitrogen due to a lack of 
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low mass stars, (2) a higher average stellar temperature in M82 than in the 
Galactic H II regions, or some combination of (1) and (2). Either mechanism 
points to the nucleus of M82 as an active site of massive star formation. 

5) Both [0 III] lines and the [N III] line have asymmetric profiles 
which are indistiguishable within the uncertainties. Our data thus show no 
evidence for large-scale variations in the N ++ /0 ++ ratio or electron density 
within the nucleus. 

We are indebted to the staff of the Kuiper Airborne Observatory for 
their excellent support and to R. Rubin and H. Dinerstein for helpful 
discussions. We also thank L. Caroff, D. Hollenbach, V. Petrosian, 

J. Simpson, M. Werner, and an anonymous referee for their comments on the 
manuscript. 
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Table 1 : Line Flux Measurements and Velocities in M82 



Rest Wavelength V LSR FWH M Flux 

Species (microns) (km/s) (km/s) (10 -1 ® W/cra 2 ) 


[0 III] 51.815 A 81 ± 25 120 ± 30 6.3 ± 0.4 

B 272 ± 25 141 ± 25 3.9 ± 0.4 

Total Flux — 10.4 ± 0.6 

Velocity Difference 191 — 

[N III] 57.330 A 126 ± 30 190 ± 20 3.0 ± 0.2 

B 314 ± 30 135 ±20 1 .3 ± 0.2 

Total Flux 4.3 ± 0.3 

Velocity Difference 188 — 


[0 III] 88.356 A 122 ± 25 150 ± 10 5.8 ± 0.2 

B 299 ± 25 150 ± 10 3.4 ± 0.3 

Total Flux — 9.2 ±0.3 

Velocity Difference 177 
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Figure Captions 


Figure 1: Observed line profile of the 51.815 pm 3 P 2 +3 Pi transition of 
[0 III] in the nucleus of M82. Data were taken 6 February 1985 with a 48" 
beam in three observations totalling 36 minutes. The horizontal axis shows 
velocity with respect to the Local Standard of Rest. Error bars represent 
± 1 standard deviation of the mean. The solid line is a two-Gaussian fit to 
the data. 

Figure 2: The 57.330 pm /2 transition of [N III] measured in 4 

observations totalling 31 minutes on the night of 4 February 1985. The 
solid line is again a two-Gaussian fit to the data. 

Figure 3: The 88.356 pm 3 P 1 ->- 3 P 0 transition of [0 III] observed for a total 
of 10 minutes on 6 February 1985. The solid line is again a two-Gaussian 
fit to the data. 

Figure 4 (bottom) : Plot of N ++ /0 ++ ionic abundance ratio divided by overall 
N/0 ratio in model H II regions as a function of effective temperature of 
exciting star, as calculated by Rubin (1985). The different curves are for 
different interstellar 0/H ratios (labelled). The top plot shows the 
dependence of N ++ /0 ++ + N/0 on the 0/H ratio for several temperatures. 
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